Shape diagram of vesicles in Poiseuille flow 
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Soft bodies flowing in a channel often exhibit parachute-like shapes usually attributed to an 
increase of hydro dynamic constraint (viscous stress and/or confinement). We show that the presence 
of a fluid membrane leads to the reverse phenomenon and build a phase diagram of shapes — which 
are classified as bullet, croissant and parachute — in channels of varying aspect ratio. Unexpectedly, 
shapes are relatively wider in the narrowest direction of the channel. We highlight the role of flow 
patterns on the membrane in this response to the asymmetry of stress distribution. 

PACS numbers: 87.16.D-,83.50.Ha 



The shape of soft bodies under flow is governed 
by strong non-linear coupling between hydrodynamic 
stresses and elastic restoring forces. The latter are of- 
ten linked with specific interface properties, like surface 
tension, or bending and shear elasticity of an elastic or 
liquid membrane. 

A widely studied flow is the confined Poiseuille flow, 
in which the behavior of red blood cells (RBCs) [iHl], 
drops [5HI], lipid vesicles [3 H EHIO], capsules [IIH16] 
or polymers [17] is often considered. Underlying motiva- 
tions include a better understanding of blood flow, or the 
possibility to manipulate these objects in microfiuidic de- 
vices for lab on chip applications. The most commonly 
reported stationary shapes are axisymmetric bullet-like 
and parachute- like shapes, the latter being characterized 
by a concave rear part. The shape of this rear part is 
very sensitive to the mechanical environment: therefore, 
observing it is a (cheap) rheology experiment in itself, as 
exemplified in Refs. [121 US], where possible membrane 
constitutive laws for capsules are discussed. Similarly, in 
Ref. [14] , the onset of curvature inversion is shown to be 
strongly dependent on the capsule's pre-stress. 

Alternatively, this shape will give indications on the 
hydrodynamic stresses on the object. It is generally ob- 
served that increasing the flow velocity, or the confine- 
ment, leads first to an increase in the fore- aft asymmetry, 
then to the apparition of a negative curvature region at 
the rear, and eventually, at least for membrane-less ob- 
jects, to break-up. Hydrodynamic interactions between 
neighboring objects are also strongly correlated to shape- 
dependent modifications of the local flow. In Ref. [4], 
small clusters of RBCs are simulated and two well sep- 
arated states (compact or loose clusters) are explicitly 
associated with two different cell shapes (shallow or deep 
parachute). This hydrodynamic cell aggregation is me- 
diated by a loop of fluid recirculation between cells, also 
called bolus [21 [3l [18] , whose apparition or disappearance 
is intimately correlated with shape changes [4 . Thus, un- 
derstanding the conditions for the apparition and stabil- 
ity of shapes and the resulting flow patterns around them 



provides a valuable entry point to build up hydrodynamic 
aggregation rules in a suspension. At larger scale, this 
self-organization of the suspension leads to specific rheo- 
logical properties like the shear-thinning of blood or the 
Farhaeus-Lindquist effect [19]. Even in the case of an 
isolated cell in a dilute suspension, shape variations will 
modify viscous dissipation, and therefore, the effective 
viscosity of the suspension. Again, specific changes in 
this viscosity can be associated with concavity changes 
at the rear and loss of membrane tension [8 . 

In this paper, we identify shape changes with flow vari- 
ations in the case of fluid vesicles, which are also sim- 
plified models for RBCs. We shall see that the general 
and therefore intuitive sketch of concavity increase with 
hydrodynamic stress must be reconsidered in the case of 
vesicles, which differ from drops or capsules by their inex- 
tensible fluid membranes. We study explicitly the effect 
of three flow characteristics: its velocity, the confinement, 
and its axial asymmetry. 

The problem is considered through experiments, nu- 
merical simulations, and theoretical calculations. In the 
experiments, we use a PDMS microfiuidic device. Vesi- 
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FIG. 1: Vesicle cross-sections in a channel of increasing aspect 
ratio, (a): v = 0.913, R = 0.35, Ca = 59. (b): v = 0.973, 
R = 0.38, Ca = 81. R and Ca are given for the square section. 
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TABLE I: Parameter ranges. For the experiments, 77 
10"^Pa.s, K = 10"^^ J and V G [10; 8300]/im.s-\ 



FIG. 2: Cross- sections for two vesicles in symmetric flow 
(exp.: R = 0.38, sim. and th.: R = 0). Left: u = 0.979, Ca = 
4191; right: u = 0.946, Ca = 3776. No fitting parameters. 



cles are prepared following the standard electroforma- 
tion method, which produces vesicles of various size and 
deflation. They are made of a DOPC lipid bilayer en- 
closing an inner solution of sugar in water. Vesicles 
are then diluted in another sugar solution and imaged 
by phase contrast microscopy. The viscosity ratio be- 
tween both solutions is close to unity. Vesicles flow 
along the x axis in a straight channel of constant thick- 
ness dz {z direction) and varying width dy {y direction). 
Their cross-section in the xy plane is observed (Fig. fl]). 
Gravity is in the x direction, so that quick centering 
is achieved [10 . Each section of given width is long 
enough for stationary shapes to be reached. The ana- 
lytical calculation is based on the decomposition of the 
shape in spherical harmonics. For simplicity, we neglect 
the influence of channel walls and use for velocity profile 

v{y^z) = V 1 — (2y/dy) — {2z/dz) , which turns out 
to be a valid approximation for low confinements. As 
in Ref. [20], solving the Stokes equations together with 
the boundary conditions at the membrane and at infinity 
and using the condition of stationary shape, we find ex- 
pressions for the amplitude of each considered harmonic 
as a function of z/, C^, and a. This method is much 
more efficient than the traditional one, which is based on 
derivation and numerical solution of the shape evolution 
equations [21 . We can now take 18 harmonics for ax- 
isymmetric case and 12 for elliptic cross-sections. The 
accuracy of the results was verified by 3D numerical sim- 
ulations using boundary integral method [22 . Finding 
theoretical and numerical approaches that fit the experi- 
mental results up to the shape details is still a challenging 
issue, as illustrated by the acute debate around phase di- 
agrams of vesicles under shear flow [22l[23]. Here, we find 
excellent agreements, as exemplified by the superimposi- 
tion of shapes in Fig. |2] 

Fluid vesicles have constant volume and surface area. 
Their deformability is directly linked to their initial 
(and constant) deflation, given by the reduced volume 
V = V/[47r(5/47r)^/V3], where V and S are the vesicle 
volume and membrane area. In the experiments, both 
are calculated thanks to axial symmetry of the shape 
in the square cross-section channel, on which we com- 
ment later on. The typical size of a vesicle is given by 
its effective radius R = y5/47r. The geometry of the 
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3 - 7 X 10^ 
10 - 6 X 10^ 



problem is characterized by two dimensionless numbers: 
the channel aspect ratio a = dy/d^ and the confine- 
ment R = 2R/ ^dydz- The capillary number, that com- 
pares viscous stress to membrane elasticity, is defined by 
Ca = V/{dydz)R^ri/ n^ where V is the unperturbed fluid 
maximum velocity, rj the fluid viscosity and n the mem- 
brane bending modulus. 

Overall, the problem is described by the four dimen- 
sionless parameters (u^a^R^Ca)- The explored ranges 
are summarized in Table [H 

Preliminary observation. — When vesicles flow from 
narrow to wide sections (the thickness dz being constant), 
their in-plane section becomes less asymmetric between 
the front and the rear, and the negative curvature re- 
gion at the rear, if any, disappears (Fig. fl]). The shapes 
are symmetric about the Ox axis [24]. A given vesicle 
flowing several times back and forth in the square chan- 
nel always shows the same 2D shape, which indicates 
that it adopts, at least, the square symmetry. Thanks 
to digital holographic microscopy, we have recently seen 
that the vesicle's transverse cross-sections are elliptical 
[26] . We then reach the a priori nontrivial conclusion 
that the vesicle shape obeys axial symmetry in a square 
geometry. We observe two types of axisymmetric shapes, 
namely parachute and bullet. We first investigate their 
existence domains in symmetric channel {a = 1) and then 
explore the influence of asymmetry for low constraints 
{R<0.5,Ca <500). 

Effect of Ca and R in symmetric channel. — For 
weakly confined vesicles {R < 0.5), clear separation be- 
tween the bullet and parachute domains is achieved in the 
(z^, Ca) space, although R varies by a factor 4 (Fig. [3[a)). 
A completely unexpected observation is the crossover 
from parachute to bullet by increasing Ca for vesicles 
with a reduced volume between 0.95 and 0.97. The same 
trend is theoretically observed for unbounded flow. The 
full 3D simulations support these results. We observe, 
however, a slight shift in Ca between experiments and 
theory. This could be attributed to uncertainties on the 
bending modulus (not measured) and to effects of ther- 
mal fluctuations and confinement, which were neglected 
in the theory. We have checked that taking into account 
a possible spontaneous curvature of the membrane has 
only a minor effect on the results. 

While, for weak confinement (i^ < 0.5), the parachute 
to bullet crossover depends on the capillary number but 
not on the confinement, the situation is completely op- 
posite for more confined vesicles. As shown in Fig. p[b). 
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FIG. 3: (color online) Shapes for ot — \. Dots refer to experimental data and full line to theory and simulations, (a): bullet- 
parachute phase diagram in the (z/, Co) plane for not too confined vesicles (0.12 < R < 0.5). Pictures in inset {u = 0.964, R — 
0.24, Ca — 20 and v — 0.959, R — 0.25, Ca — 667) illustrate the concavity change as Ca increases, (b): bullet-parachute phase 
diagram in the (zy, R) plane. For R > 0.4, Ca goes from 15 to 70000, while for R < 0.4, values are limited to the Ca < 100 
range due to shape dependency on Ca in the unbounded case. Pictures in inset (u = 0.938, R = 0.49, 0.69 and 0.89) illustrate 
the concavity change as R increases. 



for R > 0.5, well separated domains are found in the 
(z^, R) plane, despite the strong variations of Ca (more 
than three decades). This indicates that when confine- 
ment is large enough, its effect is dominant. Contrary to 
the intuition, the bullet shape is favored upon an increase 
of confinement. 

Effect of asymmetry. — We observed that, for a vesi- 
cle of given reduced volume z^, the 2D cross- sect ion, and 
in particular the concavity, is independent of R and Ca 
for low constraint {R < 0.5, Ca ^ 500), and depends 
only on asymmetry a [27 . Due to this (R^Ca) indepen- 
dence, the full 3D shape of a vesicle under flow of aspect 
ratio a is deduced from its two in-plane shapes in the 
sections of aspect ratios a and l/a, where R and Ca are 
different. For instance, in Fig. Ilia), the shapes in sec- 
tions a = 0.79 and a = 1.27 ^ 1/0.79 can be seen as the 
top and side views of the same vesicle in a rectangular 
channel of aspect ratio 0.79. We define ac as the aspect 
ratio of the channel where the crossover from concave to 
convex 2D shape occurs. If ac is larger than 1, as in Fig. 
[ija), then for 1 < a < Ofc, the in-plane cross-section of 
the vesicle is concave, and its orthogonal cross-section, 
for which we consider 1/a < olc^ is also concave. The 3D 
shape is therefore a parachute. If a > olc^ the in-plane 
shape is convex, while the orthogonal shape is concave, 
since 1/a < ac- Such a "croissant" shape has been seen 
for drops between two infinite planes [6^. If a^ is lower 
than 1, as in Fig. [lib), the same reasoning shows that 
the 3D shape is a bullet for 1 < a < l/o^c and a crois- 
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FIG. 4: Shape diagram in asymmetric channel at low con- 
straint. Exp.: each point corresponds to a vesicle with 
R < 0.5 and Ca < 500 when a = 1. Th.: R = 0,Ca< 1000. 



sant for a > 1/ac. Some weakly deflated vesicles show 
no concave shape in the explored a range, and are there- 
fore always bullets. Combining all these data, we deduce 
the general phase diagram in the (i/, a) space (Fig. [4J re- 
stricted to a > 1 for symmetry reasons) , which is in good 
agreement with the theory for i? = and low Ca- For 
Of = 1, one gets a parachute or a bullet shape; increas- 
ing the asymmetry of the channel leads to the croissant 
shape, which has a concave part in the plane of higher 
stress. 

Discussion and conclusion. — We reported the phase 
diagram of vesicle shapes in Poiseuille flow in the rele- 
vant parameter space. Our study reveals the necessity to 
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FIG. 5: Cross-sections for the same vesicle {ly — 0.985) in 
three channels, (a): dy = dz = 83/im; (b): dy = 83/im, 
dz = 92/im; (c): dy — dz — 92/im. (a) and (b) give indications 
on what would be seen in the xz plane of Fig. [lib). 
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FIG. 6: (color online) Velocity field in the vesicle co-moving 
frame. Side (a) and rear (b) view. Simulations, v — 0.95, 
Ca = 100, a = 1.2 (croissant shape). The norm of velocity is 
color-coded. Maximum velocity is about 5 times lower than 
vesicle velocity. See also [27] . 



distinguish between isotropic and anisotropic variations 
of the cross-section of the channel. While the evolution 
shown in Fig. fllis consistent with expectations (increas- 
ing confinement leads to more curved shapes), increas- 
ing isotropically the confinement leads actually to less 
curved shapes, as shown by shapes in inset of Fig. [SJ^b). 
The same conclusion is reached if the flow strength is 
increased at a given confinement. Theoretical and nu- 
merical results support these conclusions, for symmetric 
but also for asymmetric flows: as shown in Fig. [4j the less 
distorted bullet and croissant shapes are favored against 
the parachute shape upon increasing flow strength. Fig. 
[5] summarizes our main finding; isotropic expansion (a-c) 
leads unexpectedly to more shape distortion while in- 
plane-only expansion (b-c) leads to a less curved shape. 
In that case, the asymmetry between the two transverse 
directions allows to relax the constraints of perimeter and 
surface conservation at the level of a given longitudinal 
cross-section, whose evolution is therefore similar to the 
one of a drop or a capsule. Finally, the predominant role 
played by the aspect ratio of the channel is illustrated 
by the concavity increase through out-of-plane-only ex- 
pansion (a-b), which then appears to be equivalent to an 
in-plane contraction. The latter case also indicates that 
a sole 2D view in experiments might lead to biased con- 
siderations when discussing the existence of parachutes, 
asinRef. [28]. 

In asymmetric channels, vesicle shapes exhibit another 
unexpected feature: as seen in Fig. [l] and in [27 , their 
in-plane width does not increase with dy but is roughly 
constant. Therefore, contrary to drops [7J, vesicles do not 
adopt the aspect ratio of the channel but occupy a space 
in the yz plane that is delimited by a circle. The fact that 
the fluid then has to go through a small gap in the nar- 
rowest width of the channel, may lead to the intuitive but 
wrong conclusion that the vesicle should decrease its ex- 
tent in that direction (and follow more or less the channel 
aspect ratio). In order to get an insight on this peculiar 
phenomenon, we investigated the patterns of membrane 
surface flow through simulations (Fig. [g]). As soon as the 
channel is asymmetric, 4 vortices appear on the surface 



with a backward flow in the direction corresponding to 
the narrowest gap, which decreases the mean shear in the 
gap and the viscous stress on the membrane. The oppo- 
site takes place in the other direction, leading to a more 
homogeneous stress distribution. When drops adopt the 
channel cross section, all material points on their surface 
are advected backward [27 . The corresponding flow line 
possess two in-plane stagnation points at the front and at 
the back, on the main axis. This apparent singularity is 
resolved for drops through recirculation inside the drop, 
a fact that is precluded for a membrane whose material 
points must stay on the surface. This difference has al- 
ready been underlined for drops and vesicles bounded to a 
substrate and submitted to shear flow [29] . The 4- vortex 
pattern is thus the simplest acceptable flow pattern on 
the vesicle surface, if we consider that zero flow situa- 
tion is unlikely to happen for this fluid membrane under 
asymmetric constraint. This pattern can be compared to 
the one obtained recently for vesicle sedimentation [30]. 
Despite the axial symmetry of the problem, it is shown 
that a croissant shape with four vortices on its surface 
is a possible stationary shape. However, the stability of 
this solution is not discussed. 

Flow asymmetry, which can be due to the channel ge- 
ometry or to the presence of neighboring cells, will have 
therefore two consequences: from the rheological point of 
view, surface vortices of non negligible velocity (see Fig. 
[6| imply important motions of the fluid inside the vesicle 
(or the RBC), that will contribute to the net dissipation. 
These vortices, that would not be caught for instance by 
2D simulations, will also affect accordingly the flow field 
around the membrane, which will modify in return the re- 
circulating loop between cells. Whether this mechanism 
would be a stabilizing or destabilizing factor in cell clus- 
ters, that often exhibit symmetry breaking [3 , remains 
to be discussed. 
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Supplemental material for "Shape diagram of 
vesicles in Poiseuille flow" 

In this supplemental material, details of shape changes 
with channel aspect ratio a in the experiments are given. 
We also show more velocity fields on the membranes of 
vesicles and compare them with the case of droplets. 



DETAILS OF SHAPE VARIATIONS 

Due to the symmetry with respect to the x axis, we 
mainly characterize the in-plane shapes by two parame- 
ters (see Figure [t]): the concavity at the rear c = Sx/L^g,^ 
and the fore- aft asymmetry a^ = 2(5x^/Lmax- 
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FIG. 7: (color online) The different geometrical parameters 
for the in-plane cross-section of a vesicle (convex or concave). 

In the experiments, vesicles of varying reduced volume 
u and radius R were observed in channels of constant 
thickness dz and varying width dy under a flow of max- 
imum velocity V (see pictures of Figures 1 in the main 
paper). For a selection of tiny intervals of reduced vol- 
umes, we plot in Figures [sja-e) and[9]^a-d) the variations 
of as and c with a for vesicles under different confine- 
ment R and capillary number Ca that is varied either by 
changing V or because of its R dependency. In the data 
shown in Figures p[a-e) andM^a-d), R varies by a factor 
higher than 3 and Ca by a factor 100. In particular, in 
Figures [8[a) and|9[a), the same vesicle with ^ = 0.54 in 
the square cross-section channel is observed around a = 1 
for velocities V varying by a factor 100, so that Ca goes 
from 27 to 2654 in the square cross-section channel. In 
Figures ^c) andbfc), the same vesicle with R = 0.23 is 
observed on a wider interval of a for velocities V varying 
by a factor 5, so that Ca goes from 17 to 87. 

For a given v, the curves as (a) collapse reasonably well 
on a single one. In particular, as varies by a factor 5 or 
6 while a varies typically from 1 to 1.5 while the resid- 
ual variations around the master curve are much smaller 
even though R varies by a factor 3. This indicates that 
the asymmetry is controlled by the aspect ratio but not 
by the confinement (within the limit R < 0.5). In the 
range Ca ^ 500, it seems also that the Ca-dependency is 
weak. These facts suggest that the mean stress on the 
membrane (which is linked to the mean confinement or 



to fluid velocity) is not important, and that the stress 
distribution around the vesicle is the key parameter. 

Finally, figure |8|f ) shows the variations of asymmetry 
as with a for different reduced volumes. Asymmetry is 
a decreasing function of a: in-plane asymmetry is higher 
in the narrower longitudinal section of the channel. From 
the plateau observed at low a, there seems to be a max- 
imal possible asymmetry, which increases with deflation 
of vesicle, as one would expect. 

Like asymmetry, concavity is controlled by the aspect 
ratio of the channel since the curves c{a) fall down to 
at the same critical ac (concave -^ convex crossover). 
For a < OLc^ no complete collapse is observed, but this 
can be linked to the fact that tension in the concave rear 
part is weak and the membrane is therefore submitted to 
higher thermal fluctuations and might be more sensitive 
to variations in the fluid velocity V . 

A vesicle shape is also characterized by its reduced 
width t = ^xncx^/R (see Figure [7|. As shown in Figure 



10, i is more or less constant whatever a. Surprisingly, 
deflated enough {v < 0.95) vesicles are even narrower 
when the aspect ratio increases, so when the in-plane 
confinement decreases. Thus, in general, increasing the 
confinement does not lead to a in-plane shrinkage of the 
vesicle, so that the membrane is submitted to a lower flow 
stress, but only displaces the location a^ of the maximum 
thickness to the rear . 

These results imply that, at the rear, the aspect ratio 
of the transverse cross-section of the vesicle is reversed 
compared to the channel aspect ratio. It is equal to 1 
somewhere in the middle and, at the front, follows the 
same aspect ratio as the channel. These results lead to 
some comments on a previous study: in Ref. [26 we 
observed vesicles flowing in a channel of square cross- 
section. Vesicles were observed by digital holographic mi- 
croscopy, which allowed us to measure the optical thick- 
ness e{y, x) = 2An{y^ x)z{y, x), where z is the membrane 
position and An the refractive index difference between 
the inner and the outer fluid. Because a was equal to 
1, we made the assumption of axisymmetry of the vesicle 
and from this information on z we deduced values for the 
refractive index difference along the vesicle An{x). 

Strong gradients were observed, that increased with 
vesicle velocity and could lead to an increase by a fac- 
tor 2 of the refractive index difference between the rear 
and the front of the vesicle, leading to the conclusion 
of the existence of huge gradients of sucrose concentra- 
tion inside the vesicle. However, it was argued that, ac- 
cording to the permeability values found in the litera- 
ture, difference of osmotic pressure across the membrane 
should prevent such important concentration gradients, 
even though residual flow across the permeable mem- 
brane due to hydrodynamic pressure difference could ef- 
fectively lead to some sugar advection inside the vesi- 
cle (whose membrane is not permeable to sugar). The 



observed gradients in e were nevertheless attributed to 
concentration gradients because it was assumed that the 
vesicle shape should follow the channel's aspect ratio, 
even in case of small departure from the square section 
case. 

Indeed, in these experiments, poor channel quality led 
us to estimate the uncertainty on a to around 3%, that 
is, 0.97 ^ Q^ ^ 1.03. We now know that it can lead to 
important variations of the vesicle's aspect ratio along 
the X axis, because, as seen in Figure IsFe), a^ variations 
are particularly strong around a = 1, while I remains 
quasi constant. Part of the observed gradients in Ref. 
[26] could therefore be attributed to these aspect ratio 
variations. 



shown. In the axisymmetric case, no movement is seen, 
but 4 vortices appear as soon as a ^ 1, for parachutes, 
bullets or croissants. For comparison, we calculated the 
shape and surface velocities for drops with no viscosity 
contrast and capillary number Ca = 0.2. Ca is defined 
by Ca = r]VR^ /{(jde)^ where cr is the surface tension. 
Note that the choice of the value Ca value is arbitrary 
(within the condition of low enough value to prevent 
break-up), as the very different nature of the interfaces 
makes any quantitative comparison between drops and 
vesicles tricky, even with similar capillary numbers. 



FLUID MOTION AT THE SURFACE 



In Figure 11 we show additional 3D shapes for vesicles 
at Ca = 100. Tank-treading velocities at the surface are 



The flow patterns on the drop surface are completely 
different, with a front to rear movement even for a = 1, 
and no qualitative change in the asymmetric case, even 
if continuity of velocity field implies increased velocity in 
the plane of higher flow curvature. 
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FIG. 8: (color online) (a-e): in-plane asymmetry as as a function of aspect ratio a for several vesicles in the same tiny interval 
of reduced volume. R and Ca are given for the square cross-section, (a): z/ = 0.952 ± 0.004 ; (b): u = 0.965 it 0.0025 ; (c): 
ly = 0.970 ± 0.0025 ; (d): zy = 0.981 ± 0.004 ; (e): z/ = 0.987 ± 0.0025. (f): in-plane asymmetry as as a function of aspect ratio 
a for different reduced volumes (averaged over several representative vesicles with R < 0.5 and Ca < 500, from 1 to 10 vesicles 
for each curve). 
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FIG. 9: (color online) In-plane concavity c as a function of aspect ratio a for several vesicles in the same tiny interval of 
reduced volume. R and Ca are given for the square cross-section, (a): zy = 0.952 dc 0.004 ; (b): u = 0.965 ± 0.0025 ; (c): 
ly = 0.970 =b 0.0025 ; (d): z/ = 0.981 =b 0.004. The decrease of critical ac with increasing u is noticable. 
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FIG. 10: Maximum width i = ima^/R as a function of a. The initial increase can probably be associated with the proximity 
of walls (in plane confinement R/dy larger than 0.6). Note that the typical error bars on these data are around 2%, which is 
of the same order as the variations observed (but the slight decrease for high a for deflated vesicles) 
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FIG. 11: 3D theoretical shapes for vesicles and drops with fluid velocities at the surface. Color code shows the norm of the 
velocity, except for vesicles at a = 1 where there is no membrane motion and color code shows membrane tension. Vesicles : 
Ca = 100. Drops : Ca = 0.2. 



